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Abstract 



The capabilities of the functional-analytic and of the functional-integral ap- 
proach for the construction of the Hamiltonian as a self-adjoint operator on 
Hilbert space are compared in the context of non-relativistic quantum me- 
chanics. Differences are worked out by taking the one-dimensional hydrogen 
atom as an example, that is, a point mass on the Euclidean line subjected to 
the inverse-distance potential. This particular choice is made with the intent 
to clarify a long-lasting discussion about its spectral properties. In fact, for 
the four-parameter family of possible Hamiltonians the corresponding energy- 
dependent Green functions are derived in closed form. The multiplicity of 
Hamiltonians should be kept in mind when modelling certain experimental 
situations as, for instance, in quantum wires. 
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I. INTRODUCTION 



This contribution is about a basic aspect of non-relativistic quantum mechanics of a point 
particle in ci-dimensional Euchdean space M'^ subjected to a real-valued potential V{q). The 
first task in characterizing the associated quantum system is to define its Hamiltonian H 
as a self-adjoint operator acting on the Hilbert space L'^{M.'^) of complex-valued functions 
which are square-integrable with respect to the Lebesgue measure d'^q in W^. Principally, 
there are two different approaches to do so. The functional-analytic approach, see e.g. Refs. 

is concerned with constructing if as a self-adjoint extension of Schrodinger's formal 
differential operator (—1/2) A + V{q). Here A denotes the Laplacian in M.'^ and units are 
chosen such that the mass of the particle and Planck's constant h are equal to unity. In the 
functional-integral approach, see e.g. Refs. |^-^, the operator H is defined as the (negative) 
generator of a self-adjoint strongly continuous one-parameter semigroup on L^(M'^) which, 
in turn, is defined as an average of a suitable functional over the paths x : t ^ x{t) of 
a suitable Markov process in indexed by time t > 0. Under rather mild assumptions 
on the potential V, both approaches yield naturally the same unique Hamiltonian H. For 
example, this is the case if V is Kato decomposable, that is, if \V\ — V and \V\ + V are in 
the Kato and local Kato class over MJ^, respectively.^ The equivalence of both approaches is 
then concisely expressed by the standard Feynman-Kac formula 

{qt\exp{-tH}\qa)= J dW'ix] 6^''\x{t) - qt) exp I - J dt' V{x{t'))[ . (1) 

Here we have used Dirac's notation for the integral kernel of the semigroup {e~*^}t>o. 
Moreover, dW^ denotes Wiener's probability measure, which is concentrated on the set 
C(M'^, qa) of continuous paths x in starting from qa G M'^, that is, a;(0) = qa- The Dirac 
delta function accounts for the fact that the paths are conditioned to arrive at qb e R'^ at 
time t. In other words, the path average is effectively taken with respect to the Brownian 
bridge.^ 

As a prominent example we consider the inverse-distance potential V{q) = —'y/\q\, 7 G 
M. Even though the coupling constant 7 may take on both signs and this potential coincides 
only for d = 3 with the Newton-Coulomb potential on W^, we will refer to the associated 
quantum system as the (i-dimensional hydrogen atom for all c? > 1. In what follows, an 
important point is that —'y/\q\ is Kato decomposable if and only if > 2. 

The answer to the question, what happens in d = 1, will be our main concern here. In 
fact, a study of the one-dimensional hydrogen atom is of interest for at least two reasons. 
First, it nicely illustrates subtleties in the functional-analytic and in the functional-integral 
approach, which are also present in other systems with a sufficiently singular potential. 
These subtleties, while technical in nature, have important physical consequences. Only 
if the well-established mathematical techniques for handling them are generally accepted, 
one may hope for an end of a long-lasting discussion. tJiHZl Second, the one-dimensional 
hydrogen atom is of physical relevance in relation with atoms or excitons in the presence of 
strong magnetic fieldsil'tl or in quantum wires.i2Fii It is also discussed in the semiclassical 
approximation of quantum dynamics^ and as a model for an electron near the surface of 
liquid helium.ii 



2 



We will see that the standard functional-analytic approach provides a four-parameter 
family {-^b} of self-adjoint extensions of the formal differential operator (— l/2)d^/dg^ — 
7/|g|. As is suggested from the singularity of the potential, different simply correspond 
to different boundary conditions to be imposed on the wave functions ip G L^(]R) at the 
origin of the real line M. In contrast, our functional-integral approach, which seeks to keep 
the RHS of (0) finite by excluding appropriate paths, will yield only one element of the family 
{Hb}, namely Hd, the one corresponding to a Dirichlet condition at g = 0. Finally, we take 
the opportunity to rederive the early closed-form expressionll for the energy-dependent 
Green function (gb|(ifD — E)~^\qa) by a genuine functional-integral calculation. 

In order to avoid a possible confusion, we remark that the divergency of the RHS of (|1]) 
for the attractive case 7 > of the potential V{q) = ~'y/\q\ in d = 1 has nothing to do with 
a frequently encountered artifact when discretizing the functional integral inappropriately — 
sometimes referred to as "path collapse" .il This artifact occurs even in > 2. 



II. FUNCTIONAL ANALYTIC APPROACH 

In this section we present the main results of the functional-analytic approach to the 
one-dimensional hydrogen atom. The strategy we are following is the standard oneiiflil and 
was used to treat related problems or special cases. The aim is to find a dense domain 

in the Hilbert space L^(]R), on which the formal differential operator (— l/2)d^/dg^ ~ i/IqI 
is self-adjoint. 

As a first step it is often reasonable to define a formally given operator on the space 
C^(M) of arbitrarily often differentiable complex-valued functions with compact support in 
M. This is not possible here because of the l/|g [-singular ity. Therefore we start with the 
operator h := (— l/2)d^/dg^ — 'y/\q\ on the domain T>{h) := C(f (M\{0}). The operator h is 
symmetric, hence closable. The operator closure h of h has deficiency indices (2,2), as can 
be seen bv solving explicitly the corresponding eigenvalue problem in terms of Whittaker's 
functions.ca It foUowsCj that there is a four-parameter family {-^b} of self-adjoint extensions 

of h parametrizable by the self-adjoint 2 x 2-matrix B=^^^ , X, p E M. , ( E C, 

^^-"2d?"^ 

lim {27v?(g)ln(|7|g) + \ / lim \ 

<7io = B ''^ M> (3) 

lim {27(^(g)ln(-|7|g) - (^'(g)} I 1 lim (^(g) I J ' 

gTO / \ gTO / ' 

Here ^Cioc(A) denotes the set of complex-valued functions on A C M which are absolutely 
continuous on every compact subset of A, and the first and second derivatives ip' and ip" 
of if are understood to be defined Lebesgue almost everywhere. As limiting cases, we also 
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allow the entries of the matrix B to be infinite in order not to have to distinguish too many 
different cases. For example, one way to arrive at the Dirichlet condition lim^fo = 

limgjo V'(q') at the origin is by choosing ( = 0, X = p = oo. We denote the corresponding 
Hamiltonian by i^D- In the zero-coupling case 7 = 0, each describes a particle on the real 
line subjected to a certain point interaction at the origin.icl In particular, for appropriate 
B the 6- and (so-called) ^'-interaction emerge. 

Remarkably, there is a closed-form expression for the energy-dependent Green func- 
tion {qiy\ (i^B — \qa), that is, for the integral kernel of the resolvent of H^. Since this 
expression involves the Green function of Hd, we first state 



lb 



[Hb - \qa) - ^[qbqa) ^ 



X 



0(l%l - \qa\) W^,/v/=2E,i/2(v/=8^|gb|) M^/^^i/2(v/=8E|g„|) 

+ {qb < — ^ qa) ■ (4) 

Here 9 denotes Heaviside's unit-step function, V denotes Euler's gamma function and 
and W^j,^T denote Whittaker's functions.^ For g^, > 0, Eq. is just the Green function 
of the hydrogen atom on the positive half-line with a Dirichlet condition at the origin. 
Being the s-wave part of the Green function of the three-dimensional hydrogen atom (apart 
from a factor Airqbqa), it appeared for the first time in the early work of MeixnerS as a by- 
product. The second part of the next section will be devoted to a genuine functional-integral 
derivation of (H). 

For completeness and convenience of the reader we recall for the attractive case 7 > 
the twofold degenerate discrete eigenvalues and give an orthonormal basis {v^i^i , V'i^i } 
of the corresponding two-dimensional eigenspaces of Hd as obtained from the poles and 
corresponding residues of (^) 

2 

^r^ = -^, n = l,2,3,..., (5) 

^^^lil) = ©((-1)'^)^ 1^1 ^i-i(2 1^1 7/n) , A; = 1,2 . (6) 

Here Ln""* denotes a (generalized) Laguerre polynomiahlli 

Green functions of different self-adjoint extensions of a closed symmetric operator are re- 
lated via a Krein formula.ea Its particularization to the case of the one-dimensional hydrogen 
atom gives 

(g^l {Hb - E)-' \qa) = {qi>\ {Hj, - E)-' + 2^1(1 - 7/v/^))' 



x(e(g,),e(-g,))(B-.,(i5;)ll)"'(Q®f_^j)^) , (7) 
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where 

e^iE) := 27 [ln(|7|/v/^) - - -f/V^) + 2^(1)] - v^=2E (8) 
and iplz) := dlnr{z)/dz is the logarithmic derivative of Euler's gamma function. 



Since the deficiency indices of h are finite, all self-adjoint extensions possessed the 
same essential spectrum o"ess(-f^B) = [0, oo[. For the same reason and due to the fact that 
Hd is bounded below, see e.g. Eq. (j^) above, one concludes^ that they are all bounded 
below. The discrete spectrum of Hb is obtained from the poles of and consists of those 
-E G ] — 00, 0[ for which s^{E) is an eigenvalue of B 



2e,{E) = \ + p ± v/(A-p)2 + 4|CP. (9) 

A plot of £7/7 for < as a function of ^ —2E is given in Figure 1. It turns out that 
the discrete spectrum is twofold degenerate, if C = and \ = p. Not unexpectedly, for 
7 > there are infinitely many discrete eigenvalues irrespective of the choice of B. Even 
for 7 < there may be one or two negative discrete eigenvalue (s). For a given coupling 
constant 7 G M one can choose the matrix B and thus the boundary condition in such a 
way that the ground-state energy of is arbitrarily negative. 

An interesting quantity to look at is the probability current density iip{q) := 
Im{(^*(g)(y9'(g)} for (p G V{H-q). Due to self-adjointness it is continuous in q for all g G M, 
even at the origin, where it is given by 

j^(0) = Im{Clim^*(g)^(-g)}. (10) 

<jj.0 

In consequence, the choice C = implies j<^(0) = for all ip G V{H-q), which means that 
the real line is effectively split into its negative and positive half by the l/|g [-singular ity. It 
is only the resulting two-parameter subfamily that can be obtained as the orthogonal sum 
of two Hamiltonians each of which is defined separately on the Hilbert space corresponding 
to one of the half-lines. Thorough discussions of the one-parameter family of Hamiltonians 
associated with the one-dimensional hydrogen atom on the half-line have been given earlier. 



see Refs. 39-41. 



III. FUNCTIONAL INTEGRAL APPROACH 

As already mentioned in the Introduction, the Wiener integral on the RHS of for the d- 
dimensional hydrogen atom is well-defined and finite for all > 2. By now, also the explicit 
evaluationcaEl of the Laplace transform of such functional integrals is well-established. The 
one-dimensional case was explicitly considered in Ref. |5^ by formally setting c? = 1 in the 
final result of such an evaluation without any justification. The trouble with this procedure 
is that one does not know a priori to which physical problem it yields the solution. Instead, 
as a first step it is necessary to discuss the question of how the hydrogen atom can be 
reasonably defined in one dimension from the point of view of functional integration. 

To begin with, we clarify for which paths x of the one-dimensional Wiener process with 
starting point Qa E M. the integral dt' /\x{t')\ is finite. Due to continuity, one has the 
following inclusion between path sets 
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Aft{R, qa) := jx G C{R, qa) : x(t') ^ for all t' G [0, t]| 

c{.eC(E.,J:/^<oo}. (11) 



Introducing the local time0 of x at q, which is formally given by ^ [x] = dt' 6^^^ {x(t') —q) , 
one can write 



t 

dt' f dq 





Hence we deduce a further inclusion 

t 



\xit')\ J \q 



rr LtA^] ■ (12) 



ja; e C(M, ga) : ^ < oo} C ja; G C(M, qa) : Lt,o W = o} . (13) 



According to Ref. |6^, the set A/t(M, ga) equals almost surely (a.s.) the set on the RHS of 
(|13|). Therefore one infers the a.s.-equality 

t 

WY\ ^ ' ^^^^ 



which is a special case of the Theorem on p. 785 in Ref. ^ 

After this prelude the following definition for the functional integral of the one- 
dimensional hydrogen atom will suggest itself 

J dW'[x] 5^'\xit) - q,) exp | j dt' . (15) 

Five remarks are in order. 

(i) Cor. 6.3.10 in Ref. ^ implies that the restriction of the Wiener measure to A/((M, qa) 
yields again a Markovian path measure which is invariant under time reversal, namely 
the measure corresponding to the Wiener process on M killed at the origin. Using this 
fact and 



t 



lim sup [ dW^[x] [ = , (16) 

no J ^ W \x{t')\ ^ ^ 



it follows from the appropriate versions of Nakao's inequality!^ and Khas'minskii', 
lemma@ that the functional integral (|l^) exists. 
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Similar to the proof of Thm. B.1.1 in Ref. ^ the Cauchy-Schwarz inequahty together 
with (i) imply that (|15D is the integral kernel of a self-adjoint semigroup on L^(M). 
This semigroup is strongly continuous, as can be seen by adapting the proof of Thm. 
3.1 in Ref. |6^, see also Prop. 3.2 in Ref. ^ 



(iii) The generator of the semigroup defined by (|T3]) is the Dirichlet Hamiltonian Hd of 



Sec. II. This follows from the Feynman-Kac formula of Thm. 6.3.12 in Ref. and 
the fact that l/|g| is operator bounded with respect to the Dirichlet Laplacian with 
relative bound zero.0 

(iv) Restricting the paths to the set A/t(M, qa) is necessary for (|15D to exist in the attractive 
case 7 > 0. 

(v) In the repulsive case 7 < the restriction does not have any effects, if one adopts the 
convention e~°° := 0. This illustrates the fact that the existence of the unrestricted 
Wiener integral on the RHS of (|l|) does not necessarily imply that the functional- 
analytic approach leads to a unique Hamiltonian. 

In the second part of this section we will sketch how to obtain the Green function of Hd 
by computing 

00 

{q,\ {H^ - E)-' = j At e^* {q,\e~'''^ \qa) 



= j dt e^* j m^^\x] 5^^\x{t) - q,) exp \ j dt'^^ > (IT) 

for complex energies E with a real part smaller than the infimum of the spectrum of i^D- 
Thereby we follow ideas presented in Refs. ^^-p0|,|68|-[70|, where a more recent technique in 



functional integration, namely path-dependent time transformations, has been developped 
for quantum-mechanical applications. It is only with the help of this technique that a 
number of potential systems have been treated exactly in the framework of functional in- 
tegration. We are going to present only its main ideas, but all steps which are needed for 
the calculation of (0) will be explained. For a detailed and fairly pedagogical approach to 
path-dependent time transformations, the reader may consult Refs. ^, |58l - |60| , ^ . 

In what follows it will be convenient to introduce the Bessel process^ with index 1/2, 
which is a continuous Markov process taking on values in the positive half-line M+. It is 
characterized by the transition density 

/ ._ l fo^,,-i/2 / (q' - Q? \ „„„ / W + Q)' 



&r'^(g',g):=^(27rt)-i/^ (^exp|-^:^^|-exp|-^:^^|J , g', g G M+ . (18) 

Hence, Q{qq'){q/q')bl^^'^\\q'\,\q\) , q',q & M, is the (non-normalized) transition density of 
the Wiener process on M killed at the origin and the second line in (|l^) can be rewrittenlzl 
as Q{qbqa)qa/qb times 
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oo 



J dt e^* J dB,/2[x] 6^'\xit) - \q,\) exp | J dt'-^ \ . (19) 

C{K+,\qa\) lo J 

Here dBi/2 denotes the probability measure associated with the Bessel process with index 
1/2. 



The key for the calculation of (19) is the equality of measures 



d5i/2 = SidW^) . (20) 

It states that the measure of the Bessel process with index 1/2 arises as the image of the 
four-dimensional Wiener measure under the transformation S, which maps the set C(M^, r^) 
of paths y onto the set C(M+, \qa\) of paths x according to 

xit) = (Sym := [yisyit))f , rl := |g,| . (21) 

Here Sy{t) is defined to be the inverse mapping of the path-dependent time 

s 

ty{s):= [ ds'A{y{s')r, (22) 



that is, ty{sy{t)) = t. Taking (^) for granted, for the time being, we can rewrite (0) in 
terms of the paths y 

oo ( ^ \ 

(19) = A\q,\ I ds e'^^ j dW'[y] 6^^^{{y{s)f - \q,\) exp \ AE j ds' {y{s')f \ . (23) 

C(]R4,ra) I J 

In doing so, we have performed the change-of-variables t' i— >■ s' := Sy{t'), which is possible 
since ty{s) isEi dW^-a.s. strictly increasing in s, and we have used ty{oo) = oo dW^-a.s^ 
In order to calculate (|23|) we remark that 



s^'^{iyis)r '\^>'\) = -fj d'^n s^'\yis)-n) , (24) 

where J d^flrt (•) stands for the integration with respect to the usual rotational invariant 
measure on the unit sphere in over the three angles of rfc G in spherical co-ordinates. 
Using the well-known explicit expression for the integral kernel of the semigroup generated 
by the four-dimensional isotropic harmonic oscillatoioi with spring constant u;^ := —8E, 
the RHS of (p3D can be cast into the form 

oo 

72 /,,,n2 



l(^)7-w-{-i'i-i^|-|>™} 





X / d^Qr, exp \ — - n-ra \ . (25) 
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Applying the formulaEl 

[ d^ilr, exp{zn-ra} = P"!", I Ii{z\n\\ra\) (26) 
J z\rb\\ra\ 

{z G M, ly denotes the modified Bessel function of the first kind@ with index v) and per- 
forming the Laplace transformation,^ the expression ( |T9| ) is eventually seen to be identical 

to 



(2|gb|/cj|ga|) r(l - 27/u;) \ Q>{\qh\-\qa\)W2.^/uo,l/2{.^\qh\)M2^i/u,,l/2{.<^\qa\) + < ' qa) 

(27) 

Since the Green function is analytic0 on the complement of the spectrum of i^o, the desired 
result @) follows from (p7|) , provided that (^) holds. 

Following Ref. ^ the validity of the equality of measures (pOD with S as in (|21| ) is implied 
by the equality 



00 



y"dte^* j dB,/2[x\5^'\x{t)~\q,\) 



C(K+,|g„|) 



00 I 1 

4|g,| j As j m'[y] 5(1) ({y{s)f - \q,\) exp \ AE j As' {y{s')f \ , (28) 



C(R4,ra) 



which is nothing but the free-particle version 7 = of (|23|) . 

To verify ( p8D we reduce ( pT] ) for the case 7 = with the help of Ref. |78| to 

1^ (exp{v/^||g.|-|gj|}-exp{v^(|g.| + |gj)}) , (29) 

which is indeed the Laplace transform!™! of the transition density b[^^'^\\qh\, IqillL 

Within the theory of stochastic differential equations there are other ways@'li'0 to prove 



the equality of measures ([20|). In essence, they rely on Ito's formula and on a theorem due 
to Dambis, Dubins and Schwarz.0 

To summarize, the path-dependent time transformation (^Tj) relates the energy- 
dependent Green function of the one-dimensional hydrogen atom, which is the Laplace 
transform of the non-Gaussian functional integral (|T^), to the Green function of a four- 
dimensional harmonic oscillator (p3D, that is, to the Laplace transform of a Gaussian func- 
tional integral, which is explicitly known. The role of the spring constant is thereby played 
by the energy of the hydrogen atom and the role of the energy of the harmonic oscillator is 
played by the couphng constant 7. 

A path-dependent time transformation similar to (^) was used in previous functional- 
integral calculations for the radial Green function of the hydrogen atom in two or more 
dimensions. 0iiS'ii'ii Further aspects of the relation between the hydrogen atom and the 
harmonic oscillator in the framework of the stationary Schrodinger equation may be found, 
e.g., in Refs. P3,B3. 
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IV. DISCUSSION 



Clearly, mathematics alone cannot tell which particular element of the four-parameter 
family of Hamiltonians {-^b} should be chosen to model a given experimental situation. 
Instead, additional physical information is necessary. In this context it is useful to know 
that, for example, the sequence of Hamiltonians corresponding to the cut-off potentials 
— 7/(|g| + a) converges in uniform resolvent sense to the Dirichlet Hamiltonian Hd as a | 0.0 
Another fact which makes Hd stand out is that the domains of all other self-adjoint exten- 
sions contain wave functions for which the expectation values of both the kinetic and the 
potential energy are infinite — though their sum, the total energy, is always finite. Among the 
domains V{H-b) it is only T>{H-d) on which one can define the kinetic energy — (l/2)d^/dg^ 
as a self-adjoint operator. In this light, it is no surprise that the functional-integral ap- 
proach in Sec. Ill is only able to yield Hd. Indeed, there it was required to average the 
Feynman-Kac functional exp{— dt' \^(a;(t'))} with respect to a path measure generated 
by the kinetic energy. By abandoning this requirement, it is shown in Refs. p|, p3| , p4| that the 
functional-integral approach can also provide self-adjoint extensions other than the Dirich- 
let Hamiltonian. The price one has to pay is a suitable regularization of the potential near 
the singularity and a subsequent limit. For an interesting construction of a family of path 
measures corresponding to {-^b} in the zero-coupling case, see Ref. |^. 

So far, there is no systematic functional-integral analogue of the calculus of deficiency 
indices in the functional-analytic approach. To develop such an analogue would be an 
interesting programme for future research. 
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FIGURES 



FIG. 1. Plot of £7/7 for < as a function of ^jyj —2E^ where is defined in Eq. (P). The 
value of approaches 4-0(1) — 2 In 2 = —3.695 ... as 7/ \/ —2E —00. 
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